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1. Introduction 

Many important gases strongly absorb in the mid-IR (>2.5 fim) 'signature' region due to 
their fundamental rotational-vibrational transitions in this wavelength range. Optical spec- 
troscopy in the mid-IR region has potential for such applications as trace gas detection p], 
remote chemical sensing [2], and human breath analysis [3H5]. For example, human breath 
is known to contain more than 500 different 'biomarker' volatile organic compounds and 
quantification of these gases may have clinical applications. 

Simultaneous detection of several gases requires a suitable broadband or a widely tun- 
able CW source. Compared to narrow bandwidth sources, a broadband source coupled with 
Fourier transform methods offers advantages for spectroscopy including massive parallelism 
of data collection and elimination of the need for wavelength tuning. Optical frequency combs 
are particularly attractive broadband sources for spectroscopy [6], owing to their extraordi- 
nary coherence over broad bandwidth. This property has led to applications including trace 
gas detection, molecular fingerprinting and dual comb spectroscopy [7HTU] . Even though the 
resolution needed for gas spectroscopy may be much lower than the line spacing of fre- 
quency combs, one can utilize the comb structure to improve the sensitivity of trace gas 
detection. This can be done by coupling the frequency comb into a high-finesse Fabry-Perot 
cavity [TTHI3], where effective path lengths of several km can be obtained |14j . 

Synchronously pumped (sync-pumped) optical parametric oscillators represent an attrac- 
tive way of generating mid-IR frequency combs suitable for molecular spectroscopy [T^HTTj . 
Recently, our group implemented a new method for generating broadband mid-IR combs, 
based on a doubly resonant, degenerate sync-pumped OPO, which rigorously downconverts 
and augments the spectrum of its pump frequency comb [T8irT9] . Exceptionally large paramet- 
ric gain bandwidth at degeneracy combined with extensive cross mixing of comb components, 
resulted in extremely broad (> one octave) instantaneous mid-IR bandwidth extending the 
wavelength range beyond 6 fim [19]. Here we show that such a broadband source, combined 
with intracavity spectroscopy, becomes a powerful tool for trace molecular detection. 

We perform molecular spectroscopy using two such sources. One source is a periodically- 
poled lithium niobate (PPLN) based OPO, pumped at 1.56 /xm by a femtosecond Er-doped 
fiber laser [18], and the other source is an orientation-patterned GaAs based OPO, pumped 
at 2.05 jum by a Tm-doped fiber laser [19]. Both OPOs operate near degeneracy (central 
wavelength is twice that of the pump) - to obtain broad instantaneous bandwidth. In this 
work the OPO cavity itself is used as an enhancement cavity to increase the effective path 
length [20H22] . Intracavity spectroscopy of methane, carbon monoxide, formaldehyde and 
several other gases is performed by injecting gas directly into the OPO enclosure, or by using 
an intracavity gas cell with Brewster windows. We observe significant effective path length 
enhancement due to the intracavity action. In addition, we find that the measured spectral 
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line shapes may have dispersive features. Such features have previously been observed with 
cavity-enhanced frequency comb spectroscopy [TUE3], and in intracavity spectroscopy with 
sync-pumped OPOs and mode-locked lasers [T91I2HI25] . The measured spectra are compared 
to a simple model, based on the intracavity passive loss and round-trip dispersion, and 
excellent agreement between theory and measurements is found. 

The remainder of the paper is organized as follows. In Sec. [U we introduce a theory on 
how intracavity dispersion affects the spectrum and phases of the OPO comb lines. Section [3] 
presents a simple theory of intracavity spectroscopy with broadband OPOs. In Sees. H] and [5] 
we describe the experimental setup and measurement procedures. Sections |6] and FT] describe 
correspondingly results with the Er- and Tm-pumped systems. Finally, in Sec. El we discuss 
the results and make conclusions in Sec. EO 

2. Intracavity dispersion and OPO comb width 

In the case of femtosecond pumping by a mode-locked laser, the pump and the OPO output 
are represented by a manifold of comb lines sharing the same spacing, equal to the pump 
pulse repetition frequency / rep . The pump field is represented by a comb of frequencies 
v v , n = fcEO, P + nf rcp , while the OPO idler and signal follow v ifl = /ceo,* + Ure P and i/ Sjm = 
/cEO,s+ m /re P , correspondingly. Here I, m, and n are integers, and /ceo is the carrier-envelope 
offset frequency (0 < /ceo < /rep)- 

From photon energy conservation it follows that for integers I and m, there should be an 
integer n, such that 

+ V sm = Up^, (1) 

from which it follows that there are two solutions for /ceo of the signal and idler: 

/ceo,s + /ceo,? = /ceo,p (2a) 

/cEO,s + /cEO,i = /cEO.p + /rep (2b) 

For a degenerate OPO, where the signal and idler become indistinguishable, /ceo,s = /ceo,;, 
hence Eq. (j2J) becomes 

, /CEO.p / q \ 

Jceo,s — — - — , [oa) 

r /CEO.p /rep /„i \ 

/CEO,s = ^ H — • (aD) 

In order for a broadband degenerate OPO to oscillate, it is necessary that for each OPO 
comb frequency, the accumulated phase delay per roundtrip, A0 = kL, is an integer of 2-k 
(here k is the average value of the wavevector and L is the roundtrip cavity length). This can 
be achieved in vacuum (k = 2txv/c) by choosing L = c/ f rcp when /ceo,p = (c corresponds 
to the speed of light). However due to dispersion in a real cavity, A0 is no longer linear with 
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Fig. 1. For a broadband OPO to oscillate, it is necessary that for each comb 
line, Vs^m, the phase delay kL per roundtrip is an integer multiple of 2n. Due to 
dispersion in the cavity, the roundtrip phase (thick line) is no longer a straight 
line. As a result, peripheral comb lines acquire extra phase, which prevents 
them from oscillation. Dashed vertical lines correspond to the alternative set of 
OPO comb lines shifted in frequency by / rcp /2. The frequency spacing between 
the comb lines is highly exaggerated for clarity. 



v (Fig. [TJ thick line) and peripheral comb teeth acquire extra phase preventing them from 
oscillation. 

Figure|2]plots the calculated extra phase accumulated per roundtrip near the degeneracy of 
an erbium (1560 nm) laser pumped PPLN OPO, at different cavity lengths near the optimum 
sync-pumped condition. Here we took into account dispersion originating from (i) 0.5-mm 
PPLN, (ii) 2-mm ZnSe dispersion compensator, and (hi) dielectric mirror coating. These 
components have a group delay dispersion of -287 fs 2 , 296 fs 2 , and 35 fs 2 , respectively, at 
3.12 /jm. Effectively, Fig. [2] is a residue between the two curves of Fig. [TJ When L is decreased 
(from left to right in Fig. [2]), the phase curve rotates clockwise. Broadband OPO operation 
occurs when the accumulated phase for all the comb lines is an integer multiple of 2-k. 
Because there is an alternative set of OPO comb lines Eq. ( I3bl) . the interval between discrete 
cavity lengths corresponds to it instead of 2ir, in terms of the phase shift (this is expressed 
by dashed horizontal lines), so that the oscillation peaks are spaced by AL rs A pump [26J. 
Cases (a)-(c) of Fig. [2] correspond to degenerate OPO operation, with the broadest spectrum 
expected for (c), while (d) is a non-degenerate (but still doubly resonant) case, when the 
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Fig. 2. Calculated extra phase per roundtrip near OPO degeneracy (vertical 
dotted line, 96 THz, 3.12 /im) at different cavity roundtrip lengths L near the 
optimum. The length steps between (a)-(b) and (b)-(c) are 1.56 /im A pump ). 
Dispersion originating from a 0.5-mm PPLN, 2-mm ZnSe, and a dielectric 
mirror was taken into account. The cavity length is decreasing from (a) to (d). 



comb lines are clustered around separate signal and idler peaks (at ~70 and 120 THz). 

Thus, intracavity dispersion is one of the main factors that restrict the OPO bandwidth 
and can also introduce frequency chirp of the OPO pulses. 

3. Model for intracavity molecular spectroscopy 

Now we consider a high-finesse synchronously pumped femtosecond OPO in a steady-state 
condition. The round-trip evolution of the electric field A of a resonating OPO comb line at 
the frequency v S)Ta is given by 

A = tA + AA. (4) 

Here we assume that A = A Syrn exp(2iriu Syrn t), t = \t\ exp(— iA<p) is the complex round-trip 
amplitude transmission coefficient of the OPO cavity, AA is the nonlinear optical gain due 
to the presence of the pump field, and A0 - is the phase delay acquired per round-trip. The 
rationale for using the form Eq. (jlj) is as follows: For a 3-wave interaction (v s .m + = v p,n) 
in the slowly varying envelope approximation, the equation describing the evolution of the 
OPO field A(v s ^ m ) in a nonlinear crystal is given by ^A(^ Sjm ) = iKA(v p>n )A(vij)* , where 
k is a constant, proportional to the nonlinear susceptibility x^ 2 \ while assuming perfect 
phase- matching. In the sync-pumped OPO case, the nonlinear interaction results in a cross- 
coupling between large manifolds of pump and OPO comb lines (Fig. |3]), such that for the 
m-th mode 

—A{v Sim )=wJ2 My P ,n)A*{vu). (5) 
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Fig. 3. In a broadband OPO, parametric gain at a given frequency is the result 
of cross-coupling between large amounts of pump and OPO comb lines. Two 
such paths are shown. 



Here A(y Psn ) refers to the pump, A(yn) - to the 'complementary' OPO wave, and a star 
denotes its complex conjugate. Given the extremely broad spectra for both the pump and 
the OPO, the presence of narrow molecular absorption features in the spectra will be averaged 
out in Eq. (131) . Based on this reasoning, we make the approximation that the right hand side 
of Eq. §5$) is not significantly affected by weak, narrow molecular absorption features, and 
integration over the crystal length gives the constant field 'gain' AA in Eq. (Hj). 

3. A. Intracavity absorption enhancement, real amplitude transmission 

Let us first regard a simple case of high-finesse cavity when the OPO comb lines are in exact 
resonance with the cavity, exp(— iA(j>) = 1, and the round-trip amplitude transmission t(y) 
is real: t = exp(— So), where So <C 1 is the round-trip amplitude loss, such that So = S^/2, 
where <5q is the intensity (or power) loss. Also, let us assume that a small real (as in the 
center of a molecular absorption line) roundtrip absorption S\ is introduced inside the cavity, 
such that Si S , corresponding to the amplitude transmission coefficient t\ = exp(— S\). 
In this case, Eq. (J4]) results in 

AA AA _ AA 

~ 1 - t h ~ S + S ± ~ 5 (l + <5iA))' 

The relative field change SA/A, due to molecular absorption, is — Si/(S + Si) ~ — Si/S , 
instead of —S± if one measures transmission change due to molecular gas outside the OPO 
cavity. Consequently, the high-finesse OPO cavity 'amplifies' the small absorption by a factor 
1/Sq- This can also be expressed in terms of the photon lifetime inside the cavity. The photon 
number decays in a cavity as n ~ n Q exp(— 2Sof Tcp t), and thus n reduces by a factor 1/e after 
N = 1/(25q) cavity roundtrips. Hence the absorption enhancement due to the intracavity 
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action, which results in a corresponding increase of the effective path length through the gas, 
is 2N. 



3.B. Intracavity absorption enhancement, complex amplitude transmission 

When an OPO comb line acquires extra phase A0 per round trip, the transmission t is 
complex: t = exp(— 5 — iA<p). Also, molecular absorption (we assume for simplicity a 
Lorentzian line shape, which is valid at pressures > 0.1 atm), can be characterized close to 
resonance by a complex transmission function [27] 



-7 ^7-7- 



[vq -v)h- 

Here 5\ is the amplitude loss per roundtrip at the resonant frequency u , and 7 is the line half 
width. For a known gas and known path length L mo \, 81 for a given molecular transition can 
be found from the HITRAN database through the formula S\ = \8[ = Icrmoi^o^moi-^moi = 
|^n mo iL mo i, where S is the pressure- independent absorption intensity for a given energy 
transition, that is absorption cross-section per molecule integrated over the frequency [28J, 
CmoK^o) is the peak absorption cross-section, and n mo \ is molecular concentration. Similar to 
Eq. (jfJJ), we can write an expression for the field amplitude: 

AA AA 

A 



toh 1 - ^P (So ~ *A<P - 



When Acj) 7^ 0, derivative-like features appear in the spectrum, as illustrated by Fig. SJ 
Physically, dispersion in the OPO cavity introduces a mismatch A<fi(v) between the OPO 
frequency comb and the cavity resonances. Molecular dispersion imposes an additional phase 
shift near absorption resonances. This causes a change of the spectral line shapes depending 
on the mismatch. 

Accordingly, our simple theory for intracavity spectroscopy with femtosecond OPOs pre- 
dicts dispersive features in the absorption spectrum, similar to those observed with frequency 
comb spectroscopy enhanced with an external cavity [231129] . 

4. Experimental setup 

The first OPO in our experiment was based on a PPLN crystal and was pumped by an 
Er:fiber laser (Toptica, 350 mW average power, 1.56 /j,m wavelength, 80 MHz repetition 
rate, 85 fs pulse duration), producing an output centered at ~ 3.1 /im. The length of the 
PPLN crystal was 0.5 mm or 0.8 mm and a ZnSe 1-degree wedge pair was used for dispersion 
compensation and for out-coupling, as shown in Fig. |5j The second OPO was based on a 
0.5-mm-long orientation-patterned GaAs (OP-GaAs) crystal, and was pumped by a Tm:fiber 
laser (IMRA, 600 mW average power, 2.05 /im wavelength, 75 MHz rep. rate, 93 fs pulse 
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Fig. 4. When the roundtrip phase differs from an integer times 2ir, dispersive 
features of both signs appear in the molecular spectrum, (a) Calculated relative 
spectral intensity \A\ 2 near molecular resonance vs. normalized frequency at 
A(p = and A0 = 0.15. (b) Same at A0 = and A0 = -0.15. 

duration), producing an output centered at 4.1 /xm. The OPOs and their coherence properties 
are described in detail in Refs. [HI [HI [26] . Both OPOs were placed in Plexiglas enclosures 
and the pump lasers were free running. The output power of the OPOs was some tens of 
mW. 

Oscillation occurs at a discrete set of cavity lengths, separated (in effective roundtrip cavity 
length) by approximately one pump wavelength, due to the doubly resonant operation [26J. 
Figure [6] shows measured output power of the Er:fiber pumped OPO, with 0.8 mm PPLN, as 
a function of cavity length detuning. The cavity length was locked to one of these oscillation 



8 



Femtosecond 
pump laser 



Length 
stabilization 





^0 



Gas cell 



DM 



(OP-GaAs) 



Z7 



OPO enclosure 



\ 



ZnSe 
(CaF 2 ) 



FTIR 



Fig. 5. Schematic of the degenerate broadband OPO. The pump beam was 
introduced through the in-coupling dielectric mirror DM. The other five mir- 
rors are gold coated. A pair of wedges made of ZnSe (Er:fiber-pumped system) 
or CaF2 (Tm:fiber-pumped system) was used for (i) dispersion compensation 
and (ii) beam out-coupling. The nonlinear crystal was AR-coated (PPLN) or 
placed at Brewster's angle (OP-GaAs). 
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Fig. 6. Output power of the Er:fiber pumped OPO, with 0.8 mm PPLN, as a 
function of cavity length (roundtrip) detuning. 



peaks using a feedback loop including a piezo stage attached to one of the cavity mirrors. 
Figure [7] shows measured spectra for the Er:fiber OPO for a selection of oscillation peaks, 
with N 2 -purging of the OPO enclosure. The spectra are shown for the two different PPLN 
lengths used in the experiments. We observe that the maximum spectral width is limited 
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Fig. 7. Measured output spectra for the Er:fiber pumped OPO using (a) 0.8 
mm PPLN and (b) 0.5 mm PPLN. The spectra are shown for oscillation peaks 
No. 2-4 in (a) and peaks No. 1 and 4 in (b). Also shown (solid gray lines) are 
the theoretical PPLN parametric gain spectra for the two crystal lengths. 

by the parametric gain bandwidth of the nonlinear PPLN crystal. The shape of the OPO 
spectra is strongly dependent on (i) intracavity dispersion and (ii) on which oscillation peak 
the OPO is locked to [19]. Figure [8] shows measured spectrum for one of the oscillation peaks 
for Tm:fiber pumped OPO. 

5. Measurement 

We have performed spectroscopic measurements of six different molecules. Measurements of 
methane, formaldehyde, ethylene and acetylene were performed using the Er:fiber pumped 
OPO, while the absorption spectra of carbon monoxide and isotopic carbon dioxide ( 13 C02) 
were measured with the Tm:fiber pumped OPO. Measurements of the absorption spectra 
were carried out by injecting a controlled amount of gas directly into the N 2 -purged Plexiglas 
enclosure, or by using an intra-cavity gas cell with length 48 cm and a volume of ~ 30 cm 3 
(used for formaldehyde and carbon monoxide). The Brewster windows in the gas cell were 
made of 1-mm-thick ZnS that was chosen because of its comparatively low dispersion in the 
2-5 /im range (with zero-dispersion at 3.7 /im). 

The main advantages with injecting gas directly into the OPO enclosure were simplic- 
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Fig. 8. Measured spectrum for one oscillation peak for the Tm:fiber-pumped 
OPO. Also shown (gray line) is the theoretical gain spectrum of the 500 jum 
OP-GaAs crystal. 

ity and the ability to exploit the entire cavity round-trip length as a physical path length 
through the gas. The disadvantages with this approach were the large volume of the Plexi- 
glas enclosure (estimated to be 103 liter) and possible contamination of the optical surfaces 
of the OPO components. The advantages using the gas cell were the small volume of the 
cell, and the ability to easily control the flow rate and the pressure of the gas. The main 
disadvantage with the gas cell was the short cell length compared to the cavity round-trip 
distance (~ 0.5 m vs. ~ 4 m), but this could in principle be improved by using a different 
resonator geometry allowing longer intracavity gas cells. 

All spectra were measured at 1 atm. pressure and a temperature of 21.5 °C. The flow- 
rate of gas through the intra-cavity cell was 0.2 1/min. for formaldehyde, and a steady- 
state gas concentration was reached after ~ 5 min. of flowing. No flowing was used for 
carbon monoxide. Unless otherwise noted, the OPO output spectra were measured with a 
commercial (Nicolet 6700) FTIR spectrometer with a liquid N 2 -cooled HgCdTe detector. We 
used the maximum available resolution of 0.125 cm -1 of the FTIR instrument for most of the 
measurements. The total measurement time was in the range of 30 seconds up to 4 minutes, 
which includes averaging over 8 to 32 scans. First, the OPO was locked to an appropriate 
oscillation peak during measurements using the dither-and-lock technique. Then reference 
spectra were obtained by filling the OPO enclosure and gas cell, respectively, with N 2 only. 
Finally, the trace gases were injected and the absorption spectra were determined from the 
ratio between these two measurements. The free-running pump lasers resulted in drifts of 
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the OPO spectra on a time scale of the order of minutes, affecting the baseline in the relative 
OPO spectra. We therefore apply a baseline correction (determined by parts of the spectra 
between the absorption lines) of up to a few percent of the measured spectral intensity. In 
addition, an offset in wavenumber of up to ~ 1 cm" 1 was applied to the measured data, to 
compensate FTIR instrument errors and to match the positions of the measured absorption 
peaks to the HITRAN data. 

6. Results obtained with the Er-pumped system 

6. A. Methane (CH A ) 

A controlled amount of 0.2-1% methane in N 2 was injected into the N 2 -flushed OPO en- 
closure. Based on the volume of the enclosure, the concentration of methane in the OPO 
cavity was estimated to be 8.5 ppm. Figure El^a) shows measured and calculated intracavity 
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Fig. 9. (a) Measured (black) and calculated (gray) absorption spectra for 
8.5 ppm methane in nitrogen at 1 atm. pressure, corresponding to peak No. 4 
in Fig. EJ The calculated spectrum is inverted and offset for clarity, (b) Phase 
shift A<f)(v), which was used for the calculation in (a). 

spectra, corresponding to absorption of methane for the OPO using 0.8 mm of PPLN. The 
measurement time was 57 s. The round-trip phase A0 is determined from the calculated 
round-trip dispersion of the OPO cavity up to a linear frequency dependence. We adjust the 
linear part of Ac/> to optimize the fit between the measured and the calculated spectrum. 




12 



The calculated absorption of methane was based on the line intensities S(v) and line half- 
widths 7 from the HITRAN database [28]. The resulting complex susceptibility was used 
as an input to Eq. (jSJ) together with the estimated round-trip phase A0, which is shown 
in Fig. |9^b). We estimated the round-trip loss to be 25%. We observe that the simulated 
spectrum reproduces the measured spectrum well. Figure [10] shows measured and calculated 
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Fig. 10. (a) Measured (black) and calculated (gray) absorption spectra for 
8.5 ppm methane in nitrogen at 1 atm. pressure for an adjacent oscillation 
peak (peak No. 5 in Fig. [6]), compared to Fig. [9j The calculated spectrum is 
inverted and offset for clarity, (b) Phase shift A(f)(u), which was used for the 
calculation in (a). 

spectra for an adjacent (in cavity length) OPO oscillation peak (No. 5), compared to Fig. O 
In the calculation, we assumed the required cavity length change between the two oscillation 
peaks to be 1.71 /im. Nominally, the change in cavity length between two adjacent peaks is 
equal to the pump wavelength |18j . which is 1.56 /im. However, the observed spacing between 
the oscillation peaks varies within ±15%, which we attribute to intracavity dispersion. 

Based on the measured absorption spectrum and the noise of the detector and the laser 
source, we estimated a detection limit of 1.7 ppb for methane with the current experimental 
parameters. Compared to single peak detection, broadband spectroscopy detects multiple 
absorption peaks simultaneously, which leads to an increase of the sensitivity, as pointed 
out in Ref. [T7j . For example, if there are N equally strong absorption peaks, the multi- 
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line advantage improves the detection limit by a factor V N. For methane, the sensitivity is 
improved by approximately one order of magnitude, depending on the spectral resolution [TT] . 
In Sec. El we discuss the procedure of calculating the detection limit in our experiment. 
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Fig. 11. (a) Measured (black) and calculated (gray) absorption spectra for 
56 ppm methane in nitrogen at 1 atm. pressure. The calculated spectrum is 
inverted and offset for clarity, (b) Phase shift A(f)(u), which was used for the 
calculation in (a). 



Figure [TT] shows measured and calculated intracavity spectra for 56 ppm methane, with 
the Er:fiber-pumped OPO using 0.8 mm PPLN and close to zero intracavity dispersion. The 
estimated round trip loss is about 30%, which is due to increased out-coupling, compared to 
Figs. (jUD and (TTU1) . The spectrum was measured using a DTGS (deuterated triglycine sulfate) 
detector at room temperature, with 4 min. measurement time. As for Figs. (Q and (fTOj) . the 
linear part of the round-trip phase is adjusted to fit the calculated spectrum to the measured 
spectrum. 

The concentration of methane in standard air is approximately 1.8 ppm, and many of 
the methane lines are overlapping with the absorption lines of water vapor. Figure [12] shows 
a part of the OPO spectrum with no N 2 -flushing of the cavity (relative humidity 33.6%, 
temperature 21.5°C). The selected part of the spectrum corresponds to a region where a 
strong absorption line of methane is not significantly interfering with absorption lines of water 
vapor. We clearly see the absorption feature of methane in air, which is compared in the same 
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Fig. 12. Measured absorption spectrum for ambient air (black line) containing 
a methane peak. The HITRAN simulation for methane (solid gray line) and 
the three simulated (HITRAN) water vapor peaks (dashed gray line) are also 
shown . The effective path length is taken to be 8 times the round-trip length 
of the OPO cavity for the calculated spectra. 

figure to a HITRAN simulation for 1.8 ppm of methane and intracavity enhancement factor 
of 8. A reasonable match between experimental and calculated spectra is consistent with a 
round-trip loss of 25%. Simulated water vapor peaks (HITRAN) for the above atmospheric 
conditions and enhancement factor are shown as dashed lines. 

6.B. Formaldehyde (CH 2 0) 

Measurements of the absorption spectrum of formaldehyde were carried out using a gas cell 
with a continuous flow of 100 ppm formaldehyde in N 2 . In these measurements we used a 
500-/im-long PPLN crystal to obtain a sufficiently wide spectrum, corresponding to peak 
4 in Fig. [7(b), to cover the absorption region of formaldehyde from 3.3 to 3.7 /im. Figure 
ITBT a) and (b) show measured and calculated spectra of formaldehyde when the gas cell is 
placed outside (2 min. measurement time) and inside (30 s measurement time) the OPO 
cavity, respectively. The absorption was enhanced by a factor of about 6 when the gas cell 
was placed inside the OPO cavity, giving a detection limit of approximately 0.31 ppm. This 
enhancement corresponds to an intracavity round-trip loss of 33% for the signal. This loss is 
higher than for the measurements with methane, but the increased round-trip loss may be 
due to a small misalignment of the gas cell. 

For the fitted data we assumed a round-trip phase shift of A0 = 0, because this gave a 
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Fig. 13. Measured (black) and calculated absorption spectra (gray) for 100 
ppm formaldehyde in nitrogen at 1 atm. pressure. The calculated spectra are 
offset and shown on an inverted scale for clarity, (a) Extracavity spectra, (b) 
intracavity spectra. The effective path length is taken to be 6 times the length 
of the gas cell for the calculated spectrum in (b). 



reasonable overlap with the measured data. 

6.C. Simultaneous measurement of acetylene (C2H2) and methane (CH4) 

This experiment was performed to measure acetylene on one hand but also to show the 
capabilities of the intracavity OPO spectroscopy to detect more than one trace gas at once 
with reasonable high resolution. The use of a syringe allowed a precisely controlled injection 
of 400 ml of 1000 ppm acetylene in N 2 , followed by 15.0 ml of 1% methane in N 2 , into the 
N2-flushed OPO enclosure, resulting in an estimated 3.8 ppm acetylene concentration and 
1.4 ppm methane concentration, respectively. Figure H4T a) shows the measured reference 
spectrum (gray) and, underneath, the absorption intracavity spectrum (black), covering the 
absorption region of methane and acetylene, with the OPO operated with 0.8 mm of PPLN 
crystal length resonating at its peak No. 4. Figure [T4T b) indicates the dispersive absorption 
spectrum of methane and the corresponding calculated spectrum. Figure [T4T d) contains the 
detected absorption spectrum for acetylene which shows a similar amount of dispersion. Tak- 
ing both traces simultaneously did not allow optimizing the spectrum for lowest dispersion 
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Fig. 14. (a) Showing reference intracavity spectrum (gray) and, underneath, 
the intracavity spectrum (black) with the absorption features present, while 
detecting methane and acetylene simultaneously inside the OPO. (b) Exper- 
imentally measured methane spectrum (black) at a concentration of 1.4 ppm 
and the corresponding calculated spectrum (gray), (c) Round-trip phase shift 
for the calculated methane spectrum, (d) Experimentally measured acetylene 
spectrum (black) at a concentration of 3.8 ppm and the corresponding calcu- 
lated spectrum (gray), (e) Round-trip phase shift for the calculated acetylene 
spectrum. The calculated spectra in (b) and (d) are offset and shown on an 
inverted scale for clarity. 
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across the relevant wavelength ranges while maintaining a good signal to noise ratio for both 
gases. Setting the FTIR to 17 scans per trace with a resolution of 0.125 cm -1 resulted in 
an acquisition time of 60 s for the reference as well as the absorption measurements. The 
simulated spectra were obtained using the phase shifts in Fig. fl4T c) and (e) for methane 
and acetylene, respectively, while assuming a cavity round-trip loss of 25%. We estimate a 
detection limit of 0.11 ppm for acetylene with the current experimental conditions. 

6.D. Ethylene (C 2 H A ) 

A carefully controlled injection of 5.2 1 of 1000 ppm ethylene in N2 into the N2-flushed 
OPO enclosure resulted in an estimated ethylene concentration of 48 ppm along the OPO 
path. Figure [ToT a) provides the obtained absorption spectrum, resulting from the injected 
ethylene taken within 60 s, with identical settings for the FTIR as described in Sec. 16. Ci The 
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Fig. 15. (a) Experimentally measured ethylene spectrum (black) at a con- 
centration of 48 ppm and the corresponding calculated spectrum (gray). The 
calculated spectrum is inverted and offset for clarity, (b) Phase shift Ac/>(^), 
which was used for the calculation in (a). 

simulated spectrum was obtained using the phase shift shown in Fig. [ToTb) and assuming a 
cavity round-trip loss of 25%. We estimate a detection limit of 0.32 ppm for ethylene with 
our setup. 
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7. Results with the Tm-pumped system 



The Tm-system represents an attractive system because of its wide bandwidth from 2.6- 
6.1 fim at the -30 dB level. Initial measurements of intracavity spectra of water vapor and 
isotopic C0 2 in standard air have previously been presented in Refs. [T91I30] for this system. 



7. A. Carbon monoxide (CO) 

Carbon monoxide absorbs around 4.7 /iin, and the Tm-system was tuned to this wavelength 
region by selecting an appropriate oscillation peak, similar to the one shown in Fig. |HJ 
The gas cell was filled with 50 ppm CO in He and placed inside the OPO cavity. Figure 
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Fig. 16. Measured (black) and calculated (gray) absorption spectra for 50 ppm 
carbon monoxide in helium at 1 atm. pressure. The calculated spectrum is 
offset and shown on an inverted scale for clarity. The effective path length is 
taken to be 7 times the length of the gas cell for the calculated spectrum. 



[TBI shows the measured and calculated absorption spectrum of CO, with a measurement 
time of 2 min. In the calculations we assume a round-trip phase shift of A<p = 0, which 
provided a reasonable fit to the measured data. The detection limit for CO is estimated to 
be 0.27 ppm. The increased noise properties of the Tm-system, compared to the Er-system, 
are attributed to instabilities (mode hopping) in the 790-nm laser diodes used for pumping 
the Tm:fiber amplifier. We observe that the intracavity absorption is enhanced by a factor 
of ~ 7, corresponding to a round-trip loss of approximately 30%. However, this is higher 
than the expected round-trip loss of approximately 20%, estimated in a previous publication 
of the Tm-system [19], and might be due to clipping of the signal beam introduced by the 
intracavity gas cell. For a perfectly aligned gas cell, the clipping is estimated to be a few 
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percent. 



7.B. Isotopic carbon dioxide ( i3 C02) 

Isotopic ( 13 C02) carbon dioxide is naturally present in the atmosphere at concentration of 
approximately 4 ppm in volume [HI], and can be easily detected using our Tm-laser based 
OPO without purging the cavity. The resolution of the spectrometer in this experiment 
was 0.5 cm -1 (we used model 80251 Oriel-Newport mid-IR spectrometer). Figure PTT1 shows 
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Fig. 17. Measured (black) and calculated (gray) molecular spectra of isotopic 
( 13 CC>2) carbon dioxide. The simulation is based on the HITRAN database 
and is inverted for clarity. 



measured and calculated (HITRAN, with matching spectral resolution of 0.5 cm -1 ) spectra 
of isotopic carbon dioxide. To match the two spectra in their peak values, we multiplied 
the HITRAN spectrum by 15 (assuming A(p = 0), which means that our cavity absorption 
enhancement factor is larger than the one we typically observe. This discrepancy can be the 
result of an elevated (as compared to the average) amount of CO2 in the lab. 

8. Discussion 

The main motivation of using intracavity spectroscopy is to increase the effective path length 
through the gas, without the need of locking to an external Fabry- Perot cavity, or using an 
external multi-pass cell. The effective path length in our experiments dependents on the 
round-trip loss in the OPO cavity. We observe an absorption enhancement up to a factor of 
8, when the gas concentration was controlled, corresponding to a round-trip loss of 25%. This 
loss is comparatively high, and could be due to clipping of the resonating beam in the OPO 
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cavity. Based on the mirror and coating specifications, a cavity round-trip loss of 5 — 10% 
should be obtainable, resulting in an enhancement factor of 20 — 40. Maximum effective path 
lengths are thus in the range 100 — 200 m. We emphasize that once the cavity round-trip 
loss and dispersion has been determined (e.g. using a known gas with known concentration), 
both the enhancement factor and expected lineshapes are known for all molecules to be 
measured and thus absolute concentration measurements can be performed. Table [TJ shows 



Table 1. Estimated detection limits. 



Gas name 


Det. limit 


Meas. time 


Phys. path length 


System 




[ppm] 


N 


H 




Methane (CH 4 ) 


0.0017 


57 


3.75 


Er 


Formaldehyde (CH 2 0) 


0.31 


30 


0.48 


Er 


Acetylene (C 2 H 2 ) 


0.11 


60 


3.75 


Er 


Ethylene (C 2 H 4 ) 


0.32 


60 


3.75 


Er 


Carbon monoxide (CO) 


0.27 


120 


0.48 


Tm 


Isotopic carbon dioxide ( 13 C0 2 ) 


0.0024 


90 


4.00 


Tm 



estimated detection limits for each of the six molecules investigated. For calculating the 
molecular detection limit in our experiment, we have used the following procedure: First, 
we obtain the experimental set of points for the intracavity molecular absorbance S{v n ). 
Second, from the known cavity dispersion and HITRAN database, we obtain the theoretical 
absorbance spectrum T{y n ) - for the same uniformly spaced set of points v n . We remove the 
baseline from theoretical spectrum, since the baseline in our experiment is not well defined 
(we fit it to the valleys between absorption peaks). Next, we 'project' our experimental 
spectrum onto the theoretical one, that is we integrate the product S{v)T(v) over frequency, 
such that Ji = S(v n )T(v n ). Then we estimate a standard deviation a (noise) of our 
experimental spectrum and computer-generate random 'noise' N(v n ) over the set u n , with 
normal distribution and standard deviation a and project the noise onto the theoretical 
spectrum, i.e. integrate N{v)T{v) over frequency and get J 2 = iV(z/ n )T(i/ n )- The sum I 2 
fluctuates around zero with each computer run, and by doing many (> 1000) such runs, we 
get the standard deviation AJ 2 . Finally, we calculate the minimal detectable concentration as 
the concentration used in the experiment, divided by a factor ii/A/ 2 (which is an effective 
signal-to- noise factor). We think that sensitivity enhancement due to our 'matched filter' 
approach is totally consisted with the 'multiple peaks' advantage described in Ref. [17] . 
During our experiments, the pump laser was free-running, and we observed a drift of the 
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OPO spectra on a time scale of the order of minutes. This resulted in baseline drifts in 
the measured absorption spectra. One improvement would be to lock the pump laser to 
an external frequency reference. The OPO output would then be stabilized due to the self 
phase-locking mechanism of degenerate OPOs [32] . 

9. Conclusions 

By using synchronously pumped OPOs operating around degeneracy, we obtain ultra- 
broadband mid-IR radiation suitable for coherent spectroscopy in the Fourier domain. A 
large instantaneous bandwidth of up to 800 cm -1 allows detection of several trace gases si- 
multaneously. Spectroscopic detection of six molecules in trace amounts has been performed 
in the wavelength range of 2.5 — 5 /im. By injecting the gases inside the OPO cavities we 
obtained substantial enhancement of the effective path length and achieved detection limits 
down to part-per-billion level in volume. Our intracavity sensing approach offers great sim- 
plicity and compactness, which might be a great asset for future applications. The dispersive 
spectral features that we observe at some circumstances are well reproduced using a simple 
model for propagation in a dispersive Fabry-Perot cavity. These features can be predicted 
a priori from knowledge of dispersion of the intracavity elements; on the other hand, if the 
dispersion of the individual components is not known, the overall dispersion can be precisely 
mapped by injecting trace amounts of known molecules and analyzing the shapes of spec- 
tral peaks. By decreasing the OPO loss and increasing the finesse of the OPO cavity we 
expect an improvement of detection limits, down to sub-ppb levels. This will also require 
better dispersion compensation, e.g. using chirped dielectric mirrors. Utilizing intracavity 
cells with low (~0.1 atm.) gas pressure will allow better specificity of molecular recognition 
because of sharper spectral features. With further development, this system may find impor- 
tant applications in trace gas detection and real time human breath analysis - with further 
enhancement possible through the use of dual-comb multi-heterodyne methods. 
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